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Abstract The measuring stations of a geophysical network 
are often spatially distributed in an inhomogeneous man- 
ner The areal inhomogeneity can be well characterized by 
the fractal dimension Dh of the network, which is usually 
smaller than the euclidean dimension of the surface, this lat- 
ter equal to 2. The resulting dimensional deficit, (2 — D//), is 
a measure of precipitating events which cannot be detected 
by the network. The aim of the present study is to estimate 
the fractal dimension of a rain-gauge network in Tuscany 
(Central Italy) and to relate its dimension to the dimensions 
of daily rainfall events detected by a mixed satellite/radar 
methodology. We find that Dh — 1.85, while typical sum- 
mer precipitations are characterized by a dimension much 
h>. ' greater than the dimensional deficit 0.15. 

o 
en 

(^ ■ The distribution of a geophysical network is a multi-stage 

'""' . decision process which mainly relies on economic and de- 

. . ' mographic interests and on access problems in remote areas. 



1 Introduction 
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Although an ideal network of stations should be spatially ho- 
mogeneous and sufficiently dense to discriminate the mini- 
mum wave-length of the investigated geophysical phenom- 
ena, the irregularity and sparsity of observation points imply 
interpolation errors when reporting data on a regular grid. 
The areal clustering of point- sets can be measured by sta tis- 
tical indices as pointed out bv lOuchi and Uekawal (119861) or, 
when the inter-station distances are scale-invariant, it can be 



well c haracterized thanks to a fractal analysis (iMandelbrot 



1982). 



If the point-set is self-similar, i.e. any small part of it is 
the magnified version of the whole set, the set is a frac- 
tal and it can be characterized by its fractal dimension Dh, 
which is a real number with Dh < De, where De is the 
standard euclidean dimension of the embeddi ng space (in 
our case, De = 2). In Uterature several works ( Korvin et al 



199(i iLoveiov et allll986l: iMazzarella and Tranfaglial 200(A 
Olsson and Niemczvnowicj Il996t iTessier et all 119941) deal 



with the fractal characterization of a single-point observa- 
tion network and sometimes this analysis is use d as a method 

to dri v e an optimal en largement of the network (IMazzarella and Tranfaghal 
20001) . iLoveiov et all (11986) state that any sufficiently sparsely 



distributed phenomena having a fractal dimension smaller 
than the dimensional deficit 5=2 — Dh of the observing 
network cannot be detected by the network itself. Since the 
sparse precipitating phenomena are the most intense and po- 
tentially severe, they are of prominent interest, particularly 
when the network of measuring stations are constantly used 
for civil protection purposes. The aim of this study was to 
compute the fractal dimension Dh of the rain-gauge network 
belonging to the Centro Funzionale Regionale in Tuscany 
(Central Italy) and to compare Dh with the fractal dimen- 
sion d of daily rainfall events occurring in the same area, 
using an independent network for rainfall, for the month of 
July 2010. We find that Dh — 1.85, therefore giving a di- 
mensional deficit 8 ~ 0.15. On the other hand, all rain pat- 



terns give a fractal dimension d > 0.6, well above 5. A rough 
extrapolation of data for li as a function of 24-hours rain 
thresholds suggests that our rain-gauge networks might fail 
to record precipitation events whose intensity is about 75 
mm/day or more. 

The paper is organized as follows. In section |2] we de- 
tail the data used in this study and the methodology adopted 
to evaluate Dh- In section [3] the computation of the fractal 
dimension of the rain-gauge network is presented and com- 
pared to those obtained for rainy days. Finally in section |4] 
results are discussed with reference to the potentiality and 
limits of the applied methodology. 



2 Data sources and methods 

2.1 Data sources 

The location of the rain-gauges belonging to the Centro Fun- 
zionale Regionale (CFR) in Tuscany (Central Italy) is shown 
in Figure[T] Its establishment has been a long-term decision 
process involving several local institutions over more than 
20 years. The network comprises 377 stations and encloses 
several basins over an area of about 23000 km^ (yielding a 
density of about one station every 60 km'). More than 90% 
of the stations are located below 810 meters. The biggest 
inter-station distance (in other words the size or diameter of 
the point-set) is about 250 km. These data r riake the geogra- 



consist of a mosaic of instantaneous surface rainfall intensi- 
ties (SRI) with a spatial and temporal resolution of 1 km and 
15 minutes, respectively. Several daily precipitation amounts 
were tested; thresholds values of 1 mm, 2 mm, 5 mm, 10 
mm, 15 mm, 20 mm, 25 mm, 30 mm, 35 mm per day were 
used to analyze the different phenomenology linked to pre- 
cipitation, from weak to moderate regimes. 



2.2 Methodology 

While euclidean geometry deals with ideal geometric forms 
and assigns dimension to points, 1 to lines and so on, frac- 
tal geometry deals with non-integer dimensions. The frac- 
tal, or Hausdorff, dimension Dh has been the most common 
used measure of the strangeness of attractors o f dissipative 



dynam ical systems that exhibit chaotic behavior ( iGrassberger and Procacci 



1983bl) . Since for experimental data the value of D h is diffi 



cult to determine using the box-counting algorithm dStrogatz 



19941) . we computed the fra ctal dimension D2 of point -set 



using t he method proposed in lOrassberger and Procaccia ( 1983at 
I983bh. as al s o found in the literature ( Korvin et alTl 1990l). 



Loveiov et all(ll986l) . lMazzarella and TranfaghalOOOOh . lOlsson and Niemc: 
(119961) 1 In the case under examination we choo se to use D? 



as a goo d approximation of Dh, since as stated in lGrassberger and Procaccic 
(Il983bl) . D2 < Dh and inequalities are rather tight in most 



cases. 



In the present study we compute the correlation dimension 

in a 2i dimensional space but in general, to obtain D2 given 

phx^our network similar to that studied by|Mazzarella and Tranf^ljatset {X,}^; with X, € R", we have to consider the 

(| ZOOU|) . correlation integral C{R) that counts the number of pairs 

{X;,Xy} such that || X,— Xy || is smaller than a given thresh- 
old R > 0, with nil being the standard euclidean distance 
in R". In formulas: 



Satellite imagery acquired by the Meteosat Second Gen- 
eration (MSG-2) satellite in the infrared (IR) channel cen- 
tered at 10.8 jim was used in this study as a proxy to detect 
and monitor cold clouds systems. The study period is July 
2010 while the spatial and temporal data resolutions are 4.5 
X 4.5 km^ and 15 minutes, respectively. A brightness tem- 
perature (Tb) threshold was used to identify cold cloud sys- 
tems th at are most likely to be a s sociat ed with convective 
activity. iKolios and FeidasI (120071 l2010h used a Tb of 228 
K to best identify convective systems in the Mediterranean 
area based on a set of lighting data. The same tempera ture 
threshold of 228 K was used bv lMorel and Senesil (120021) for 
assessing the climatology of the Europe an MCSs, this value 
being very close to that of 22 1 K used bv lGarcia-Herera et al 
(I2OO5I) for Spain. This low-temperature threshold allows to 
investigate mostly anvil regions and embedded areas of ac- 
tive deep convection (Johnson et al. 1990). In this study, a 
228 K Tb threshold value was chosen for identifying very 
deep convective events over Tuscany. Furthermore, discrim- 
ination between precipitating and non-precipitating cold cloud 
systems, previously subjected to Tg threshold test, was per- 
formed using RADAR data provided by the DPCN (Na- 
tional Civil Protection Department) radar network. The data 
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(1) 



where is the Heaviside function and where j^ij^_ j ^ is the 

normalization factor so that C{R) tends to 1 for R tending to 

infinite. 

If the rain-gauge network is a fractal then C{R) grows like a 

power: 



C{R)o:R"\ 



that is 



log{C{R))cxD2log{R). 



(2) 



(3) 



Therefore, one can derive D2 from the regression coefficient 
of relationship (O. 

In order to determine the correlation dimension Do of the 
rain-gauge network described previously, we computed the 
correlation function defined in equation ([T]), as described in 
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Fig. 1 Location of rain gauges (right side of the picture)- belonging to the Centro Funzionale Regionale network in Tuscany (Central Italy, on left 
side of the picture). 



Loveiov et all (119861) . i.e. we determined the cumulative fre- 
quency distribution of the inter-station distances for the to- 
tal number of 377 stations. The distances were determined 
by spherical trigonometry, using geographic coordinates and 
ignoring elevations owing to the smallness of the elevation 
with respect to the two horizontal dimensions. 
F or what concerns the value s of th e parameter R, as done 
in iMazzarella and Tranfaglial (l2000h . we started in comput- 



Distribution of minimum distances 

Average 4.285 km 



ing the inter-station distances defined in equation ([T]i from 
1 km. This value was gradually increased by a factor of 1 . 1 
up to 250 km since, as expected by definition of C{R) given 
in equation ([1]), for all R > size (area of interest) the correla- 
tion integral C{R) saturates to 1 and log{C{R)) saturates to 
0. 

For experimental data the linear behavior of log{C{R)) on 
log{R) is Umited to a scaling region Sr, i. e. only for R be - 
longing to the interval Sr = [RmimRmax] (IStrogatzlll994l) . 
This happens because C{R) is underestimated from those 
points near the edge of the set so that the criteria to deter- 
mine the bounds of Sr need to be analyzed in each singu- 
lar casejLiebovitch a nd TothI 1989h. According to the liter- 
ature (IForrest and Witten Il979t iGrassberger and Procaccia 
1983al: iKorvin et all 1 19901) . the upper limit R„uix is chosen 
equal to one third of the diameter of the area (about 80 km). 
In order to choose the lower limit /?„„„, we didn't perform 
any statistical significance computation, since in our case the 
correlation coefficients are statistically significant at 99% 
confident level for all /? > 1 km. Rather, for each station 
we computed the distance of the nearest neighbor and took 
the average of this distribution as the meaningful index of 
the points separation. 

In Figure 12 we plot this distribution; the average of nearest 
neighbor's distances is about 4.2 km and this value is consid- 
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Fig. 2 Distribution of nearest neighbor's distances for each station in 
the point-set. Average value is 4.2 km which is taken as R„,in, the lower 
limit for the regression of log(C(S)) on log(S). 



ered as the lower limit /?„„■„ meaningful for the regression. 



3 Results 

The linear fitting between log{C{R)) and \og{R) within the 
scaling region Sr bounded by /?„„„ = 4.2 km and R„,ax = 
80.2 km yields a slope, and thus a correlation dimension 
value D2 of 1.85. Figure |3j shows the results of this regres- 




Fig. 3 Log-log plot of correlation integral C{R) on R with scaling 
region delimited by S„„„ = 4.2 km and R,„ax = 80.2 km. The corre- 
sponding slope of the regression line which determines the correlation 
dimension D2 is equal to 1.85. 



sion. The dimensional deficit 5 of the network, defined as 
the difference between the dimension of the embedded space 
andD2, is (2-£)2) = 0.15. 

The value of 5 should be related to the dimension d of 
rainfall phenomena. From a climatological point of view, the 
area of interest (central Italy) is mainly affected by convec- 
tive storms or frontal systems, depending on the seasonality. 
Convective storms are of uppermost interest for our analysis 
since they are smaller, more or less separate rainfall areas 
displaying a considerable spatial variability and thus suit- 
able for fractal analysis. They are typical of the warm sea- 
son (from June to September roughly). The frontal storms 
are characterized by continuous rainfall areas of large spa- 
tial extensions and are typical in autumn and winter seasons. 
Using remote sensing and ground instruments described in 
section |2] we collect data for every day in July 2010. First 
step is to select all the pixel in the MSG-2 15-minutes dataset 
having a brightness temperature below 228 K so that we 
can obtai n a point-set (i.e. pixel-set) of potential precipita - 
tion cells (IKolios and FeidasluOlOtlMorel and Senesill2002b . 
Secondly, to assign a rain amount to the selected pixels we 
use the radar data. For each selected pixel in the MSG 15- 
minutes dataset we retrieve the surface rainfall intensities 
(SRI) as estimated by the RADAR data provided by the 
DPCN (National Civil Protection Department). Finally for 
each day of July 2010 we add all the 96 daily images (for 
each day we have one image every 15 minutes) and obtain a 
daily estimate of precipitation amount. Rain estimates were 
processed in order to compute the correlation dimension us- 
ing the method detailed in section |2] In figure |4] we plot the 
correlation dimensions of rainfall events registered in the 




exp model: Y= 1.66 exp(-X/ 31.13 ) 
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Fig. 4 Correlation dimensions of rainfall events (average values and 
standard deviations) registered in Tuscany in the month of July 2010 
for different thresholds. Horizontal gray line is 5, the dimensional 
deficit of rain-gauge network. 



month of July 2010 in Tuscany. Daily rainfall events were 
divided on the basis of prescribed thresholds, chosen equal 
to 1 mm, 2 mm, 5 mm, 10 mm, 15 mm, 20 mm, 30 mm, 
and 35 mm. For each threshold. Figure |4] shows the average 
and standard deviation values of correlation dimension of 
the rainy pixel-set for those days having a significant num- 
ber of points that registered an amount of precipitation above 
the threshold. 



4 Discussion 

The present study achieved the issue to estimate the areal 
sparseness of the monitoring rain-gauge network belonging 
to the CFR owned by Tuscany Administration by means of 
the fractal (correlation) dimension D2. In Table [1] we com- 
pare this value with dimension D2 found in other, similar 
studies in the literature. Except the cases of Australia and 
Canada, where the dominance of inhabited areas along the 
coast lowers the value of the fractal dimension, our D2 value 
is in good agreement with the others. However, we have 
to point out that the computed correlation dimension D2 
must be han dled with care bec ause, according to the Tso- 
nis criterion dTsonis et alll 19941) . the minimum number N,„in 
of points required to produce a correlation integral with no 
more than an enorErr (normally Err = Q.QSDe) is approx- 



Reference 



# of points Area coverage D2 



Loveiov et al (1986^ 



— Korvin et al (1990) 

— Tessier et al (1994') 

— Olsson 

and Niemczynowicz (1996) 

— Mazzarella 

and Tranfaglia (2000) 

— Present study 



9563 

3593 

414 

-65000 

7983 

230 

215 

300 

377 



global land 
France 
Canada 
Australia 
global land 
~ 10000 km^ 



1.75 
-1.8 
-1.5 
1.42 
1.79 
-2 



: 38000 km^ 1.84 

: 38000 km^ 1.89 
: 23000 km^ 1.85 



Table 1 Comparison of D2 computed in the present study with other 
correlation dimensions found in literature. It is also reported the bib- 
liography references, the number of points taken into account and the 
extent of geographical area. 



imately 

which, in our case, means Ni„i„ — 600 whereas we have 377 
stations. 

In order to analyze the dimensional deficit 5 = 0.15, 
we need to compare it with the fractal dimension of rain- 
fall events, as done in the previous section. As for the pre- 
cipitations fallen in July 2010, even the more intense ones 
were characterized by a fractal dimension d > 0.6, much 
greater than the dimensional deficit 5 = 0.15. We stress that 
the used rain data (see Figure |4|i are independent from the 
rain-gauge network used to determine the fractal dimension 
Dh = 1.85. Therefore, these preliminary results allow us to 
state with a good confidence that our rain-gauge network 
was precise enough to record all precipitation events oc- 
curred in July 2010. Empirically we can suppose that the 
fractal dimension goes to zero as the threshold increases; 
this because intense precipitation events (at least thermo- 
convective ones) are more scattered. On the other hand light 
rains are more homogeneus and then associated to a (de- 
creasing) linear trend for small thresholds. We then choose 
to fit the values plotted in Figure|4]with an exponential model 
which is almost linear for small values of x and decreases to 
zero as x tends to infinity. We can then suppose that 



y = 



:Ae 



-X/XQ 



where the independent variable x on the x-axis represents the 
daily amount of rain, y represents the fractal dimensions and 
the parameters are A = 1.657 and xq = 31.133. The model 
is calibrated using the just the first values of x (thresholds 
from 1 mm/day up to 15 mm/day, continuos line in Figure 
lU, since above we don't have any significant statistics (just 
two days registering at least 20 pixels with a precipitation 
above 20 mm and one day registering at least 20 pixels with 



a precipitation above 25 mm). 

The exponential regression intersects 5 = 0.15 for 



xs = -J^oln ( - 



that is xo.15 ~ 75 mm/day. In other words our data suggests 
that rainfalls with daily amount equal or above 75 mm/day 
might correspond to a fractal dimension d < 8,so that these 
events could not be detected. This value is based, by con- 
struction, on the remote sensed data and ground instruments 
and on the phenomenology of rainfall events. Ongoing ef- 
forts are directed toward the improvements of the accuracy 
of instruments and toward the calibration of the algorithms. 
Further studies are required to investigate the relation- 
ship between correlation dimension D2 of the observing net- 
work (but we are rather interested in the dimensional deficit 
5) and the dimensions d of rainfall events. Firstly the most 
important improvement of the research is to expand the statis- 
tics of precipitating events considering several months for, at 
least, a couple of years. Moreover it would be interesting to 
take into consideration the fall/winter precipitations, which 
are mainly associated with cold and warm fronts, to evaluate 
the different behavior of dimension d and check if, for some 
thresholds, it drops below 5. 
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